Interaction effects and phase relaxation in disordered systems 



I.L. Aleiner 1,2 , B.L. Altshuler 3 ' 4 , and M.E. Gershenson 5 
1 Ruhr-Universitat-Bochum, Bochum, D-44^80, Germany 
2 Department of Physics and Astronomy, SUNY at Stony Brook, Stony Brook, NY 11794 
3 NEC Research Institute, 4 Independence Way, Princeton, NJ 08540 
A Physics Department, Princeton University, Princeton, NJ 08544 
5 Serin Physics Laboratory, Rutgers University, Piscataway, NJ 08854-8019 

(February 1, 2008) 

This paper is intended to demonstrate that there is no need to revise the existing theory of 
the transport properties of disordered conductors in the so-called weak localization regime. In 
particular, we demonstrate explicitly that recent attempts to justify theoretically that the dephasing 
rate (extracted from the magnetoresistance) remains finite at zero temperature are based on the 
profoundly incorrect calculation. This demonstration is based on a straightforward evaluation of 
the effect of the electron-electron interaction on the weak localization correction to the conductivity 
of disordered metals. Using well-controlled perturbation theory with the inverse conductance g 
as the small parameter, we show that this effect consists of two contributions. First contribution 
comes from the processes with energy transfer smaller than the temperature. This contribution is 
responsible for setting the energy scale for the magnetoresistance. The second contribution originates 
from the virtual processes with energy transfer larger than the temperature. It is shown that the 
latter processes have nothing to do with the dephasing, but rather manifest the second order (in 
1/g) correction to the conductance. This correction is calculated for the first time. The paper 
also contains a brief review of the existing experiments on the dephasing of electrons in disordered 
conductors and an extended qualitative discussion of the quantum corrections to the conductivity 
and to the density of electronic states in the weak localization regime. 

PACS numbers: 73.20.Fz, 73.20.Jc, 73.23.-b, 



I. INTRODUCTION 

At low temperatures, the classical conductivity of disordered conductors (normal metals and semiconductors) 
is determined by the scattering of electrons off a quenched disorder (e.g., impurities and defects). This residual 
conductivity is given by the Drude formula: 
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(1.1) 



where vf andpi? are the Fermi velocity and Fermi momentum of the electrons, correspondingly, m and n denote their 
mass and density, r is the transport elastic mean free time, while I = vft is the mean free path of electrons. 

This expression can be justified provided the elastic mean free path I is sufficiently large, since all the corrections 
to the classical Drude conductivity ( [lTl| ) are of the order of h/{ppl) <C 1. However, in low dimensions, these quantum 
corrections to the conductivity (QCC) diverge when temperature T decreases. Eventually, they drive the system to 
the insulating regime. 

QCC are of a substantial importance even for conductors that are far from the strong localization regime: in a wide 



range of parameters QCC, though being smaller than the classical conductivity (1.1), determine all the temperature 
and magnetic field dependences of the conductivity. The systematic study of QCC was started almost two decades 
ago. The comprehensive review of the status of the problem from both theoretical and experimental viewpoints can 
be found in several papersErQ. 

According to their physical origin QCC, can be divided into two distinct groups. The correction of the first type, 
known as the weak localization (WL) correction, is caused by the quantum interference effect on the diffusive motion 
of a single electron. For low dimensional (d = 1,2) infinite systems, WL QCC diverges, and this divergence should 
be regularized either by magnetic field or by some other dephasing mechanism. 

The second type of QCC, usually referred to as the interaction effects, are absent in the one-particle approximation; 
they are entirely due to interactions between electrons. These corrections can be interpreted as the scattering of 
an electron off the inhomogeneous distribution of the density of the rest of the electrons. One can attribute this 
inhomogeneous distribution to Friedel oscillations produced by each impurity. The role of the electron-electron 
interactions in this type of QCC is to produce a static self-consistent (and temperature-dependent) potential. Such a 
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potential does not lead to any real transitions between single-electron quantum states. Therefore, it does not break 
the T-invariance of the system and neither affects nor regularizes the WL corrections. We will elaborate more on this 
point in Sec. |l[ 

However, the interactions between electrons are by no means irrelevant to WL QCC. Indeed, these interactions 
cause phase relaxation of the single-electron states, and thus result in the cut-off the divergences in WL corrections. 
This dephasing requires real inelastic collisions between the electrons. When the temperature T and magnetic field H 
are low enough, this dephasing provides the leading mechanism to regularize QCC of the weak localization type. This 
means that in this regime interactions between electrons can not be taken into account perturbatively: the higher is 
the term of the perturbation theory expansion the stronger it diverges. For the reasons, which will become clear below, 
in this paper we restrict ourselves by the lowest order of the perturbation theory in the interactions. This is possible, 
provided the dephasing and, consequently, the WL effects are mostly determined by the external magnetic field. 
Under these conditions interactions between electrons modify the dephasing rate slightly, and thus cause temperature 
dependence of the magnetoresistance. As a result, we end up with two classes of interaction contributions to the 
conductivity: corrections to WL QCC due to the interactions (class I) and the genuine interaction QCC (class II). 

Recently, there were several attempts to revise this conventional picture. Based on their own measurements, as well 
as on experimental results of other groups, authors offl'Q argued that apparent saturation of the dephasing rate at 
T — > has an intrinsic nature - scattering of the electrons by "zero-point fluctuations of electromagnetic environment" . 
This suggestion was claimed to be put onto more "theoretical grounds" by Golubev and Zaikinlffl. 

Even though these suggestions (as we strongly believe and aim to explain below) can be rejected using purely 
physical arguments, they still generate a sympathy in some part of the physical community. The situation becomes 
worrisome, since the paperfl, at the first glance, looks like a respectable calculation an experimentalist could rely on. 
Since it is much easier to observe the saturation rather than not to observe it (a possible reason for this saturation is 
discussed, e.g., in Ref. |l^), one can expect that in the future, quite a number of experimental and non-experimental 
papers confirming the role of the zero-point motion in the dephasing mechanism will appear. Having all this in mind, 
we have undertaken in this paper a straightforward analysis of the effect of electron-electron interactions on the weak 
localization. We explicitly demonstrated that the dephasing (i.e. the characteristic scale of the magnetoresistance) 
is determined by the real processes with the energy transfer smaller or of the order of temperature. The role of 
the processes with the transmitted energy larger than temperature (which was referred to as "zero-point motion" 
in Refs. is to produce type two QCC. The explicit formula for those corrections in the second order of the 

perturbation theory in the inverse dimensionless conductance has been obtained for the first time. 

The paper is organized as follows. Section |h| reviews qualitatively physics of the quantum corrections to the 
conductivity. In the same section we introduce the notion of the dephasing time t v . Section III briefly outlines the 
essential content of the "new theories" tm which predict zero-temperature dephasing. SectionTTy] contains the main 
idea of our calculation and the results, which clearly contradict to all the conclusions of Refs. ||,[|~Technical details of 
this calculation are described in section |y| Section [v| is devoted to a critical analysis of Refs. g,^. We discuss their 
results and possible sources of mistakes made in these papers. In the same section we show that the results of these 
references disagree with the available experimental data in two- and three- dimensional systems by several orders 
of magnitude and the agreement in one-dimensional case claimed in Ref. ^ is a pure coincidence in the region of 
parameters were the weak localization theory is already not applicable. Our findings are summarized in Conclusion. 



II. QUANTUM CORRECTIONS TO THE CONDUCTIVITY: A QUALITATIVE PICTURE 

This section contains only rather old results and ideas. We review them here in spite of the fact that there is 
a number of review articles that discuss ia .details both qualitative and technical parts of the theory of disordered 
conductors in the weak localization regimetTB. The main reason to do this is to make this paper self-contained, so that 
any reasonably educated physicist can read it without consulting the reviews continuously, and to underline the main 
physical concepts which are being rejected by Refs. ^-|^. We also propose somewhat new intuition about the physical 
origin of the effects of the electron-electron interaction. Nevertheless, experts in the field can skip this section. 



A. Weak localization 



Anomalous magnetoresistance in disordered conductors (doped semiconductors and metals) has been recognized 
for almost 50 yearsE3. For a long time this phenomenon has remained a puzzle. The theoretical understanding of 
the anomalous magnetoresistance emerged as a spinoff of the theory of Anderson localization. It turned out that the 
correction to the conductivity, which is due to quantum interference at large length scales, is very sensitive to weak 
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magnetic fields. The quantum correction itself may be much smaller than the classical conductivity. Nevertheless, the 
weak field magnetoresistance is dominated by this correction, and its basic features (its amplitude, dependence on both 
magnitude and direction of the magnetic field, etc.) are very different from that of the classical magnetoresistance. 
Since the quantum correction can eventually drive the system to the Anderson insulator, the regime when this 
correction is small, is called the weak localization (WL) regime, and the theory of the anomalous magnetoresistance 
is now a part of the theory of weak localization. 

A qualitative physical interpretation of WL is usually based on the following arguments, see e.g. Ref. [l3|. Consider 
an electron diffusing in a good conductor, ppl 3> h. The probability w for the electron to reach, say, point i starting 
from point /, see Fig. [l], can be obtained by, first, finding the semiclassical amplitudes A a for different paths connecting 
the points, and, second, by squaring the modulus of their sum: 



5> 



0- 



(2.1) 
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The first term in Eq. (2.1) is nothing but the sum of the classical probabilities of the different paths, and the second 
term is the quantum mechanical interference of the different amplitudes. For generic pairs a, 0, the product A a A*a 
oscillates on the scale of the order of Xf — fi/pF as the function of the position of point /. This is because the lengths 
of paths a and (3 are substantially different. Since all measurable quantities are averaged on the scale much larger 
than Xf, such oscillating contributions can be neglected. There are pairs of paths, however, which are coherent. An 
example of such paths is shown in Fig. ^b. The paths 1 and 2 coincide almost everywhere but the loop segment BEB 
(see Fig. |I)d) is traversed by trajectories 1 and 2 in the opposite directions. In the absence of the magnetic field and 
spin-orbit interactions, the phases of the trajectories 1 and 2 are equal. Therefore, the contribution of this paths to 
the probability w becomes 
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\A 2 



2ReAxA* 2 = 4 \Ai 



(2.2) 



i.e. it is twice as large as the classical probability. Thus, in order to determine the value of the weak localization 
correction to the conductivity SawL, one has to determine the classical probability to find such a self-intersecting 
trajectory. 

In order to find this probability, we label all the trajectories by the time t it takes for a classical particle to go 
around the loop. The classical probability dP that the diffusing particle returns into the phase volume dV at a given 
time is 
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(2.3) 



where D = vpr/d is the diffusion constant, and d is the number of dimensions. The relevant phase volume can be 

estimated as vpdt (SpS<p) d ~ 1 , where Sp and 6<j> characterize the transverse distance between the paths 1 and 2 at the 
intersection point, see Fig. |l]b, and 6<j) is the intersection angle between them. For the interference between paths 
1 an d 2 to be effective, the uncertainty relation pp5(j)5p ~ fi should hold. We substitute this estimate for dV into 
Eq. (2.3) and sum up over all the trajectories: 
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(2.4) 



where the negative sign is due to the fact that the returning trajectory should arrive at the interse ctio n point with 
the momentum almost opposite to the initial one. At low dimensions d = 1,2, the integral in Eq. (2.4) diverges at 
the upper limit. Below we will return to the cut-off of this integral. 



The explicit calculation first performed in Ref. |14| gives the result similar to Eq. (2.4) up to a numerical factor: 

Sa W L = r / dtC(r,r;t), (2.5) 



where a is defined in Eq. (1.1) and v is the density of states per one spin. The retarded classical propagator 
Cooperon is the solution of the equation 
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DV\ )C = 8{r 1 -r 2 )5(t). 



(2.6) 
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FIG. 1. Examples of the classical (a) non-coherent and (b) coherent paths between points i and /. The region of intersection 
of trajectories is blown up in the inset. 



It follows from the previous discussion, that the weak localization correction originates from the quantum interfer- 
ence within the pairs of time reversed paths. Therefore, this correction is extremely sensitive to any violation of the 
time reversal invariance, e.g. magnetic impurities, external magnetic field, external microwave radiation etc. (we will 
not consider the spin-orbit interaction in this paper). 

Let us discuss first the effect of the magnetic field. The role of the inelastic processes will be considered in the 
followiag subsection. The effect of the magnetic field on the conductance is a manifestation of the Aharonov - Bohm 
effcctlia. Let us again consider a pair of the coherent trajectories 1 and 2 from Fig. |l|b. Since those trajectories traverse 
the loop in the opposite directions, the Aharonov-Bohm phases acquired by them are of the same magnitude, but 
opposite signs: 



Ax -> Ax exp (i-^-HS^ ; A 2 -> A 2 exp (-i-^HS 



(2.7) 



where S is the directed area of the loop and H is the external ma gnetic field. As the result, Eq. (2.2) should be 
modified as ReAi^ = |^4i| 2 cos (2\HS\ and, accordingly, Eq. (2.4) now reads 
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(2.8) 



where Ah = y/hc/eH is the magnetic length. The last transformation in Eq. (2.8) relies on the fact that the area S, 



swept by the path is the normally distributed random quantity for a diffusive system. For the infinite three-dimensional 
systems and two-dimensional films in the. magnetic field perpendicular to the film plane, (S 2 (t)} ~ (Dt) 2 . In a Id case 
(wire) with the transverse dimension aJlJ (S 2 (t)) ~ (Dt) a 2 ; the same expression holds also for a two-dimensional film 
in the magnetic field parallel to the film plane with a being the film thickness. p-. 
As the result, the weak localization correction is cut-off at the large time limit by th equal tota 
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d=3, d = 2 (n\\H); 
d=l, d = 2 (n±H), 



(2.9a) 
(2.9b) 



where n is the vector normal to the plane of the film. 

More rigorously, the weak localization correction to the conductivity is still given by Eq. (2.5), however, the equation 
for the Cooperon (|]^) is modified: 
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C = 6( ri -r 2 )S(t), 



(2.10) 



where the vector potential A due to the magnetic field describes the Aharonov-Bohm phase we have already discussed. 
For d = 3 and d — 2 (n || ff), this equation is equivalent to the one for a charged particle in the magnetic field, and 
£Ih from Eq. (2.9a) has a meaning of the cyclotron frequency. For one-dimensional systems and films (n _L H), it 
reduces to 
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C = S(n-r 2 )S{t), 



(2.11) 



with th being defined by Eq. ( |2.9b| ). 

Let us now define the dephasing time t v and the corresponding length scale L v = Dt v in a pure phenomenological 
way. First, we repeat that the phase itself is not a directly observable quantity, therefore one always has to specify the 
observable which is sensitive to the phase relaxation. As we have already discussed, the WL correction is cut by the 
magnetic field. On the other hand, it can also be cut by any process that breaks the phase coherence. We introduce 
the function Fd(x) so that 
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Function Fd(x) is different for different dimensionalities d, and it has the asymptotic behavior Fd(x) = const, i«l 
and Fd(x) oc x d l 2 ~ x , x 1. Therefore, t v is defined through the scale of the magnetic fields at which magnetoresis- 
tance occurs. Since both the conductivity 5a and time th can be measured with a great accuracy, one can extract 
t v by fitting the magnetoresistance. In order to define this time with the numerical coefficient, the explicit form of 
functions Fd(x) is needed, which requires some theoretical analysisEZl. We postpone further discussion of this point 
until the end of the following subsection. 



B. Interaction effects 



It was already mentioned in the introduction that effects of the interaction between the electrons on the conductivity 
can be subdivided into two groups. The corrections of the first group are associated with scattering of a given electron 
on the static Friedel - like oscillations of the charge density caused by disorder. These corrections have nothing to 
do with the energy transfer between the electrons, i.e., with inelastic processes. The second group originates entirely 
from the inelastic scattering of the electrons. Since the physics involved is completely different for those two groups, 
we consider them separately. 



1. Correction to the conductivity and tunneling density of states 

First, we will discuss the Hartree type corrections to the tunneling density of states and to the conductivity assuming 
that the interaction between electrons is determined by a short range two-particle potential V . Later we will show 
how to generalize obtained formulas for the Fock correction and for the Coulomb potential. The rigorous derivation 
of the results is given in Ref. § Here, our goal is to describe the qualitative physical picture. 

Consider the local density of states of a non-interacting system 

i/fcrJ^hMtOf^e-eO, (2.13) 

i 

where ipi and q are the eigenfunctions and eigenenergies of the disordered systems respectively. Hereafter, all the 
energies are measured from the Fermi energy. 

Semiclassically, the density of states ( 2.1 3| ) can be presented in the fori 



v(e, r) m const + Im^ A Q (e, r), (2-14) 

a 

where A a (e,r) is the quantum mechanical amplitude corresponding to the path a starting at and returning to the 
point r, see Fig. ||a, and "const" stands instead of a function of energy that is smooth on the scale of the Fermi 
energy. 
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The second term in Eq. (2.14) oscillates on the scale of the order of Xp — fo/pF as the function of the position of the 
point r, because the lengths of the paths a,/3,7 are much larger than A^. As the result, this contribution vanishes 
upon disorder averaging and there is no interference contribution to the density of states of non-interacting electrons 
in the disordered metals. 

Let us demonstrate now that the interactions between electrons change the situation drastically. In the Hartree 
approximation, the only effect of the interactions is reduced to the effective one-electron potential Vjj, which should 
be considered in addition to the random potential of impurities 

V H (r) = J dr' p(r')V (r - r') . (2.15) 

Here p(r) is the electron density and V (r — r') is the potential of the interaction between two electrons located at the 
points r and r'. We assume that this potential is weak and short-range on the scale of the order of the e lastic mean 



free-path I. In a clean infinite system, the electron density is homogeneous, and the Hartree correction (2.15) leads 
only to a uniform (and thus not observable) shift of the chemical potential. In disordered or finite systems, however, 
the electron density has an inhomogeneous term 

p(r) = fde n F (e) Im^A a (e;r) (2.16) 

J a 

which oscillates on the scale of the order of Xf', for a single impurity these oscillations correspond to the well-known 
Friedel oscillations. Here, rti?(e) = 1/(1 + e t/>T ) is the Fermi-Dirac distribution function. 

Scattering of an electron by these density oscillations results in the singularity in the average tunneling density of 
states. To demonstrate this, consider the scattering off the Hartree potential in the first-order perturbation theory. 
According to the usual rules of quantum mechanics, the correction to the return amplitude 8 A is given by 

6A(e; r) a J dr' V H (r') £ A (r, r'; e) A, (r', r; e) , (2.17) 

0,1 

where A a (r,r';e) is the quantum mechanical amplitude for the electron to get from point r to point r' using the 
classical path a, see Fig. ||b; A a (e; r) = A a (r, r; e). 




a, 




FIG. 2. (a) The returning path that contributes to the density of states v and the oscillating part of the electron density 
p; (b) Classical paths responsible for the correction to the amplitude of the return; (c) Three paths producing non-vanishing 
correction to the density of states. The regions of the intersections of the trajectories are blown up in the insets. 
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Substituting Eqs. (|2.15D and ( |2.16D into Eq. ( |2.17| ), we obtain with the help of Eq. ( |2.14| ) 

<fc/(e.r)w /" dr'dr"^ (r' - r") J de'n F {e') [ImA a (e', r')] [W^ (r, r"; e) A 7 (r", r; e) 



(2.18) 



Contribution of a typical set of paths (a, f3, 7) to the sum ( 2.18 ) oscillates as the function of r. However, similarly to 
the case of weak localization, there exist special sets of coherent paths, with non- vanishing contributions to the sum. 
An example of such sets is shown in Fig. |2|c. In this case, paths 7 and/3 are chosen to be parts of the trajectory a. 
We will denote them a < and a>. Let us neglect in the sum in Eq. (2. IS) all the terms except the contributions of 
the sets (a, ,3,7) of this type: {(3 — a<, 7 = a> The phase of the quantum amplitude A a of a path a is connected 
with the classical action S a which corresponds to this path: A a cx e lSa /%. Note that difference between the energies 
e and e' is small, and that points r' and r" are close due to the short range nature of the potential V. Expansion of 
S a in coordinate and energy with the use of identities 



dS a (r 1 ,r 2 ;e) 
de 
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dS a (T-i,r 2 ;e) 
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dS a (T-i,r 2 ;e) 
dr 2 
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(t a is the time it takes for an electron with energy e to go along the path a from the initial point r 2 to the final point 
while p\i\\ are respectively its final and initial momenta) allows to rewrite A a (e',r') as 



A a (e',r')=A a (e,r")e 
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As a result, Eq. fl2.18j ) acquires the form 
5v{e.r) « J2 J dr ' V ( r ') ex P i tr ' (Pi ~ Pa) ] Re / dr " J dt ' n F ( £ ') e u ^-^ A a (e, r")A* a< (r, r"; e) A* a> (r", r; e) . 



(2.19) 



To further simplify equation ( [2.19 ) , we notice that for the diffusive system, directions of the final and initial momenta 
are not correlated. Therefore, the first factor in Eq. ( 2.19] ) equals to 



J drV (r) exp [ir (p{ - p 4 Q )] -> V (p x - p 2 ), 



(2.20) 



where V is the Fourier transform of the interaction potential and the bar denotes the averaging over the Fermi 
surface. Also, since the aver aged density of states is translationally invariant, one can replace the integration over the 
interaction point r" in Eq. (2.19) by the integration over the observation point r. We use the identity 



dvA* a< (r, r"; e) A* a> (r", r; e) = tt a A* a (r"; e) , 

which can be derived using the exact relation for the Green function. The fact that the result is proportional to the 
length of the trajectory t a can be easily understood since the observation point r can be inserted anywhere along the 
path a. 



After the integration over e', we obtain from Eq. ( 2. IE ) 



Sv(e)*V( Pl -p 2 )R C J2- 



irTt a 



sinh nTt n 
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(2.21) 



As a result, the correction to the density of states is expressed in familiar terms of the return probability, which 
we have already dealt with in a previous subsection! Substituting P(t) ~ l/(Dt) d ' 2 for the return probability, we 
estimate the correction to the tunneling density of states as 



V (Pi - Pa) 



dt cos (et) — 



irTt 
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(2.22) 



Therefore, the one-particle density of states v(e) gets significantly modified due to the scattering of a gi ven el ectron off 
the charge density oscillations of the rest of the electron gas. At low dimensions, the correction Su(e) (2.22) diverges 
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at large times when T = and e = 0. One can interpret the divergence as the Bragg scattering of the incident electron 
off Friedel oscillations. Finite electron energy e sets the upper limit for the integration because the scattering remains 
resonant only as long as the length of trajectory vpi times the difference Tivp/e between the momentum becomes of 
the order of unity, which corresponds to the condition ei < h. A finite temperature T also cuts the correction: it 
smears out the oscillating contributions related to the traje ctorie s longer than hvp/T. 

One can obtain a more familiar and convenient forma of ( 2.22 ) by making Fourier transform in time 



,H 



(e) 



= Im / — 



2^ 



tanh 



2T 
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2T 
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2V 
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DQ 2 ) 
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(2.23) 



Here we introduced the short-hand notation V for the constant defined in Eq. ( 2.20| ). 

As far as the correction to the conductivity for a short-range interaction is concerned, it is directl y rel ated to the 
correction to the density of states. For instance, one can determine Hartree QCC directly from Eq. Q2.23 ) 
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(2.24) 



Inclusion of the Fock part of the self-consistent potential is also straightforward. It amou nts to the addition of the 
Fourier transform of the interaction potential — V(q = 0) to the potential 2V in Eqs. (2.23) and (2.24). The negative 
sign reflects the fermionic statistics of electrons. 

For the Coulomb interaction between electrons V(q — 0) diverges. This indicates that the screening should be taken 
into account. As a result, the correction to the tunneling density of states acquires the form 
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(2.25) 



where the dynamically screened interaction potential is given by 



U(Q,u) 



{v(Q) 



2v- 



DQ 2 



-iuj + DQ 2 



(2.26) 



For a long range interaction, there is no straightforward relation between the density of states and the conductivity 
as Eq. ( [2.24 ), since the large part of such interaction should be gauged out. One finds (see Ref. |l| and Sec. |v| of the 
present paper for the further technical details) 
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For the further references, we show here the results of integration in Eq. (2.27): 

= 2 f h 



Sa c (T) 



2zV, V 



d= 2; 



6a c (T) 



irh V 2ttT \ 2 J 1 



(2.27) 

(2.28a) 
(2.28b) 



where ((x) is the Riemann zeta-function, £(3/2) = 2.612 . . .. 

It is crucial ( and we will keep emphasizing it) that the frequency uj in Eqs. (2.25), transfered by the Coulomb 
interaction and (2.27) is larger than the temperature T. We saw that physical meaning of this correction is the elastic 
scattering of the electron from the potential (Hartree or Fock) created by the rest of electrons. These processes do 
not cause any phase breaking since they do not violate the time reversal symmetry. 



2. Energy and phase relaxation 

We turn now to discussion of the second type of effects due to the electron-electron interaction - the inelastic 
scattering. 

In early papers (see e.g. Ref. |l^) on the theory of localization, the dephasing rate 1/t v was considered to be of the 
same order as the inelastic collision rate in clean conductors. The latter can be expressed as the sum of the electron 
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- phonon (l/r e _ p ^ ~ T 3 /0|,) and electron - electron (l/r e _ e ~ T 2 /Ep) contributions, where Ep and Qd are the 
Fermi and Debye energies correspondingly (here and almost everywhere below we put fi = 1 and = 1)- Even under 
this assumption, the e — e contribution dominates at low enough temperatures. JtJpecame clear later that the static 
disorder enhances strongly the e — e contribution to the inelastic scattering rateEHlEl, while 1 /r e - p /i is less affectedcil. 
As the result, both dephasing and energy relaxation rates at low temperatures are governed by collisions between 
electrons. 

To recall the main results on the e — e dephasing rate, let us start with a single electron excitation, assuming that 
T = 0, and the rest of the electron gas occupies states below the Fermi level. The dependence of l/r e _ e on the 
excitation jenjamy e (the energy of an electron counted from the Fermi level) can be determined in a perturbative 
calculationEMB The result [Eq.(4.4) of Ref. [| can be rewritten through the dimensionless conductance g(L) (the 
conductance measured in units of e 2 /h ~ 1/25. 8kCl) of a d - dimensional cube of the size L 

1 



( - C d ——, L e = y/D/e, (2.29) 

T e-e(e) g[L e ) 

where Cd is the dimcnsion-and-coupling-constant-dependent coefficient. For a weakly interacting Id electron gas 
Cd = V%- Equation ( 2.29] ) can be also rewritten as 



—L- = C d 6 1 (L e ), (2.30) 
r e - e (e) 

where Si(L) = (i d ^) _1 is the one-particle mean level spacing in a d - dimensional cube of the size L and v is the 
one-particle density of states. P 

There are several interpretations of this result. One of themtl is based on the concept of the interaction time which 
becomes much longer in the disordered case due to diffusive rather tha n ball istic motion of electrons. It is also possible 
to appeal to statistical properties of exact one-electron wave functionsBEJ, and we outline this interpretation below. 

The inelastic rate l/r e _ e (e) is determined by a pair of collisions between electrons with all four energies - two initial 
(e a and e 7 ) and two final (e$ and eg), all the energies here are counted from the Fermi level. 

Given the typical absolute value M a 7 s of the matrix element for such a collision in a sample of the size L, the 
inelastic rate can be estimatednJ'Ell with the help of the Fermi Golden Rule 

"7 ~, 7 m3 6 I M aPlS\ 2 8 ( £ a ~ e /3 + e 7 ~ e &) [1 ~ n F{tp)] (^7) [1 ~ n F (es)] ■ ( 2 -31) 

The matrix elements can be represented as integrals of products of four exact one-particle wave functions. In 
a disordered system, these wave functions oscillate randomly in space, and are only weakly correlated with each 
other. As a result, the matrix elements are random and for L smaller than L e (Od case), their typical absolute value 
M(L) turns out to be of the order of 5\{L)/g(L), where the small factor g~ x reflects the weakness of the correlation 
between the wave functions. (In the limit g — ► oo Random Matrix Theory is valid; according to this theory, there 
is no correlation a t all between different eigenvectors and the non-diagonal matrix elements vanish.) As a result, 



integration in Eq. ( 2.31 ) for Od case can be easily performedE.il and for T = 

7m- (2 - 32) 

Thi s rat e incr eases with e, and at L — L e , it becomes of the order of i5i(L c ). This estimate corresponds exactly to 



Eqs. (2.29) and (2.3C), and it remains valid even for large samples, L > L e , since l/r e „ e cannot depend on L in this 
limit. 

This observation enables us to estimate the dependence of the matrix elements in a d - dim ensio nal sa mple on the 
energy transfer uo = e a — ep, when u> is larger than the Thouless energy. Comparison of Eqs. (2.3C) and (2.32) yields 
for this dependence 

2 S 1 (L) 3 6 1 (L U ) h{Vfl/ 2+dj2 
\M a p 1 $\ — 5 = TTd — x w ' CJ =l e a _e /3|- (2.33) 

This energy dependence of the matrix elements reflects the properties of noninteracting disordered system and is not 
sensitive to the energy distribution of the electrons. 

It follows from Eq. ( p. 33 ) that the matrix element s div erge when d < 4 and u — > 0. At T — this divergence is 



not dangerous because of the two integrations in Eq. (2.31). However, the situation changes when the temperature is 
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finite. Assuming that the energy transfer uj is much sma ller t han the temperature T, we perform the integration over 
e 7 and eg in Eq. ( 2.31 ) and obtain with the help of Eq. ( 2.33| ) 



duiH(oj); n(w) 



T 



(2.34) 



The integral over uj diverges in the infrared limit for d = 1,2. Therefore, l/r e _ e is ill-defined at finite temperatures 
and in low dimensionsc^l. 

This is not a catastrophe, though: l/r e _ e itself has no physical meaning. When the energy relaxation rate l/r e , 
(i.e. the inverse time of thermalization of an excitation with energy e much larger than temperature T) is considered, 
the quasielastic processes are not important. To be more precise, for the energy diffusion the relevant quantity is 



dujuj 2 lZ(uj) , 



which is perfectly convergent at the lower limit. Therefore Eqs. ( |2.29 ) and (2.30) give a good estimate of l/r e 

e 1 




FIG. 3. (a) Example of the coherent paths involving excitation of the electron hole pair at point ti. The phase difference 
between paths 1 and 2 is <5</>i2 = u)(t — 2ti), where uj is the energy of the excited pair. 



As it was mentioned in the end of Subsection II A, the characteristic scale of the magnetic fields in the weak 
localization is determined by the dephasing rate l/r^. It turns out that this dephasing rate at low dimensions is 
governed by the quasielastic processes. 

Let us illustrate the microscopic mechanism of the dephasing. Consider a pair of coherent paths 1 and 2 from 
Fig. ^b. Suppose now that an electron during its motion along path 1 emits (or absorbs) an electron-hole pair with 
the energy uj, see Fig. ||. Clearly, the interference term will be lost, unless the electron going through path 2 emits 
(absorbs) the same electron- hole pair. In order to evaluate, e.g., the conductivity, one has to sum up over all the 
final states. The contribution of a final state with the extra electron-hole pair contains the interference part from 
trajectories 2 and 1 with the emission of the very same pair. The only difference between paths 1 and 2 is the following: 
for path 1 the electron has momentum ki = kp + e/vp during the time ij, see Fig. [|, and then it has momentum 
fcj — uj jvp during the remaining time t — ti\ the electron going along path 2 has momentum fc, during the time t — ti 
and momentum ki — uj/vp during time rj. As a result, the phases acquired by the trajectories are no longer the same: 



Sfi2 = J kidl — J k2<il = uj(t 



2*0- 



(2.35) 



Since the effect of the dephasing is significant only for Sip 12 — 1, the processes with the frequency, smaller than the 
inverse length of the path do not lead to the d epha sing, and thus the relevant quantity for the path with length t is 
not the total electron-electron relaxation rate ( 2.34 ) but rather 



T V (t) 



duj1Z{uj) 



T 



1/' 



>Dt 



(2.36) 



It is this quantity that determines the dependence of the weak localization correction on temperature. If the magnetic 
field is strong enough, t h < t v (t h ), see Eqs. (2.S), the dephasing introduces a small correction of the order of 
th /t v (th)- On the other hand, if the magnetic field is weak, the lengt h of the relevant trajectories should be 
determined self-consistently by putting t — t v in the right-hand side of Eq. ( [2.36| ) 
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T 



g (i/Dt^ 



(2.37) 



Solution of Eq. yield J 



2/3 



T R u e 2 l ( irh 
h 2irh 11 \e 2 Rn 



d = l, 



d = 2, 



(2.38a) 
(2.38b) 



where u\ is the conductivity per unit length of a one-dimensional system (wire) and Ru is the sheet resistance of the 
two-dimensional film. 



An expression for 1/t v with the exact numerical coefficient (as Eqs. (2.38)) rather than an order of magnitude 



estimate is meaningful only if it is supplied with the explicit dependencies of the conductivity on 1/t v and H. In 
order to evaluate these dependencies, one notices that the dephasing is dominated by the processes with the energy 
transfers much smaller than temperature. In this case, instead of considering the interaction between electrons, one 
can think about an electron moving in a classical electric field, created by the rest of the electron gas. Since this field 
fluctuates in time and space, the calculation of observable quantities involves averaging over these fluctuations. The 
program was r ealiz ed in Ref. where it was shown that the weak localization correction is still given by expression 
similar to Eq. (2.5) 



2(7 f 00 

Sa WL = / dt (C (r, r; t. - t)). 



where the Cooperon in the external vector potential A is given by expression similar to Eqs. (J2.6|) and (2.1C) 



ch 



A(n,-t/2) 



C ( ri ,r 2 ;t,t') = 5 (n - r 2 ) S(t - t'), 



(2.39) 



(2.40) 



where A(r, t) is the fluctuating vector potential created due to the emission and adsorption of the real electron-hole 
pairs. Those fields are Gaussian ones with the propagator defined as 



(A i A j ) 



2T c 2 Q i Q j 

2— o— Im U{Q,uj) 



2T c 2 Q'Q 1 

<7d LO 2 Q 2 



(2.41) 



where U(Q,ui) is the screened interaction potential defined in Eq. (2.25), and ad is the conductivity in d dimensions. 
The last expansion in this formula is possible since the relevant frequencies and wavevectors are such that Q ~ 
y/bj/D ~ l/L v <C l/a s , where a s is the screening radius. Formula ( p. 41 ) is consistent with Eq. (4.21) of Ref. |l]. The 
first term in the latter formula describes the transverse fluctuation of the electromagnetic field which is usually small 
and we will not take them into account. Calculation described in great details in Sec. 4.2 of Ref. |l] yields 



5o WL = 



In 



8(TWL 



T n H \2 



1 



111 



2ir 2 Tl T (TH + T0) ' 



SfTwL = 



i 



irh [— In Ai (t v / tjj)] 



d = 2, (n|| H) - 
d = 2, (n±H); 

d = l, 



(2.42a) 
(2.42b) 
(2.42c) 



where the dephasing time t v and the time th associated with magnetic field are given by Eqs. (^38|) and @ 
respe ctivel y for different dimensions d, Ai(a;) is th e Air y-function, and ^(x) is the digamma-function. Forms of 
Eqs. ( 2.42 ) are consistent with the general relation ( 2.12 



It is important to emphasize that derivation of Eqs. (2.39) - (2.41) is neither an assumption nor a phenomenological 
recipe - it is the result of the parametrically justified procedure. All the corrections to these formulas are small as 
\/{tJT). As such, this calculation is already self-contained, however some people find it useful to interpret the results 
appealing to Gedanken experiments similar to that of R ef. |28| . 

For further consider ation s it is useful to rewrite Eq. ( 2.42| ) in a slightly different form. First, we perform a gauge 
transformation of Eq. (2.40): 
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2e 

-*Vi - -rA(n) 

cn 



+ i( X (r u t/2)-x(ri,-t/2)) 



C(r 1 ,r 2 -t,t') = 8{r 1 -r 2 )8(t-t'), 



(2.43) 



where A{r) represents external magnetic field, while the fluctuating electric field cause d by the interaction between 
electrons enters through a fluctuating scalar potential x(r,i). According to Eq. (2.41), correlation function of this 
scalar potential in (cj, Q) representation equals to 



(x 2 



2T , N 2Te 2 
lmU(Q,oj ) = — FS . 



(2.44) 



Secondly, for the regime th t v , the fluctuating potential % can be taken into account by perturbation theory. This 
yields the correction to the conductivity of the form 



SawL = 
C{Q,u) = 



a 



d d Q 
(2n) d 



C(Q,0) 



d d Qd d q f T dcu fATe 2 



(2tt 



,2,/ 



_ T 2tt \ a d q 2 



-ito - 



DQ 2 + ^~ 



([C(Q, 0)] 2 C(Q + q,u)- \C(Q, u;)\ 2 C(Q + q, 0)) 

(2.45) 



where the first term in brackets is just a weak-localization correction for non-interacting electrons and the second 
term is the perturbative effect of dephasing on the weak localization. Since in d = 1 the integral over u) is determined 
by its lower limit, the upper limit \u>\ = T is not important. In d = 2 the integral is logarithmical and the exact 



definition of the c ut-off is also not necessary. Integration in Eq. ( 2.45 ) gives exactly the same result as the direct 
expansion of Eqs. (2.42) in the small parameter tjj/t^. Here, and through the rest of this paper we will consider only 



one-dimensional wires and the films in the parallel magnetic field to avoid unnecessary technical complications. 



C. Summary of this section. 

This section is a review of the results on the effects of electron-electron interaction on the conductivity of disordered 
metals. We tried to give the physical interpretation to all of the rigorous formulas. The main message is that the 
microscopic mechanism of dephasing is the excitation or absorption of real electron-hole pairs. That is why the 
frequency transferred through the interaction propagator in the formulas for dephasing can not be much larger than 
temperature, since excitations with higher energy are forbidden by Pauli principle. As a result the dephasing rate 



I/t^ vanishes when the temperature tends to zero (see Eqs. ( 2.3q )). The interaction correction associated with higher 
frequency are not due to the real processes but rather due to the scattering of the electron off the static self-consistent 
potential created by all the other electrons. Therefore, this scattering has nothing to do with the dephasing. In the 
next section, we outline "new theories" of dephasing which attempt to revise this physical picture. 



III. RECENT "THEORIES" OF ZERO-TEMPERATURE DEPHASING 

Recently, Mohanty and WebbEl suggested the following expression [formula (4) of Ref. [?J for the dephasing time 
1 e 2 f , f d d k f duj ^ f %w\ 

T, 



, du I t / — Tkjcoth — fc 2 exp (— Dk \u\ — iuju) . (3.1) 

ah 2 J J (27T) d J 2tt \2T) K ' 



We will denote this time t p to distinguish it from the results of previous section. We also use integration £J^ d instead 
of their . ( We do not quite understand why the authors of Ref. |t] conclude that "In a quasi-lD system, for 

any two interfering paths the wave vector k must be a two-dimensional vector"). Formula (3J) was obtained as an 
extension of the derivation of Refs. ^7],[l] in the form popularized in Ref. [50| This extension can be hardly called a 
derivation. _ _ 

Golubev and ZaikinQ'B claimed that they derived Eq. fl3.l| ). To obtain Eq. (fO| ), one needs just to substitute Eq.(74) 
of Ref. H into Eq. (71) of the same paper, see also Eqs. (7) and (8) of Ref. g. We will discuss derivation of Ref. g in 
more details later. 

Beforehand we would like to to emphasize the differences between Eq. (3.1) and the results of previous section. 
Though authors do not write it explicitly, the magnetic field dependence of the conductivity in Refs. is presumably 
given by 
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5a 



WL 



d d Q 



1 

TH 



(3.2) 



Consider the strong magnetic field th <C t p limit. As it was already mentioned, in this limit one can expand Eq. (3.2) 
in small parameter th/t p . Let us restrict ourselves to the first order of this expansion, substitute r p from Eq. (3.1) 
and integrate over u in Eq. (3T). The result is 



Sa 



WL 



a 



d d Q 



(27T) 



- d C{Q,Q) 



d d Qd d q [ 1/T dto [2Te 2 \ / hoj\ t „ fn tl2 



-4W ■ 



(3.3) 



where C(Q,w) is defined in Eq. 

Let us now compare Eq. (3.3) with 



(E-4ED . 

the conventional result ( [2.45 ). The first terms of both equations describe 
weak localization of non-interacting correction, and they simply coincide. There is no dramatic difference in the 



contributions to the second terms from small energy transfers. Indeed, the integrand in the second term of Eq. (3.3) 
at %ui <C T is of the same order of magnitude as the one in the second term of Eq. (2.45). In order to avoid an 
unphysical conclusion that some dephasing can be caused by the fluctuations with ui = or q — 0, i.e., by a static or 
uniform potential, the authors of Refs. ||-^| introduced an ad hoc infrared cut-off at ui ~ 1/tj, 



The difference at high frequencies is much more serious. While Eq. (2.45) limits the frequency integration over the 
energies smaller than temperature, Eq. ( |3.3| ) extends the integration over the whole range of frequencies. As a result 
of this extension, the dephasing rate becomes finite at T — > 0. Thus, according to Eq. ( |3.3| ) the processes with the 
energy transfer much larger than temperature do break the time reversal symmetry and do cause dephasing. On the 
other hand, all the physical arguments of Sec. || apparently contradict to this statement. 

We emphasize that it is extremely difficult to accept the idea of the zero-point motion dephasing from the qualitative 
point of view. Consider a quantum particle which moves in the environment of the harmonic oscillators: each oscillator 
is characterized by its frequency oj. The result of the collision of the particle with an oscillator depends on the relation 
between T%u> and the temperature. If hio <C T, the oscillator is in an excited state with a high probability. The particle 
energy, counted from the Fermi level, is also of the order of T, i.e. much larger than huo. As the result, the energy 
conservation does not prevent the energy transfer between the particle and the oscillator, i.e. the probability for the 
inelastic collision is substantial. 

The situation for fuj 3> T is quite different. Indeed, up to the exponentially small terms, the oscillator is in the 
ground state, and the particle has the energy smaller than hoj. Therefore, the energy transfer is forbidden by the 
energy conservation and the collision is elastic. Therefore, there is no difference whatsoever between the collision with 
such an oscillator and with the quenched disorder, which definitely does not dephase. As to the zero-point motion, 
their energy is simply the difference between the ground state energy of the oscillator and the bottom of the harmonic 
potential. Since, this is the energy of the ground state, it can not be absorbed, and the smearing of the wave-function 
of the oscillator is just a renormalization of the cross-section for the elastic processes. 

However, maybe our intuition is wrong? We did not find any convincing physical argument which would require to 
reconsider the common understanding of dephasing in disordered systems in Refs. |6|-|9| 

The only theoretical argument in favor of the zero-temperature dephasing presented up to now is the calculation 
of Ref. ||. We do not want to reject these calculations entirely on the basis of the qualitative picture. The history 
of Physics knows examples when a calculation preceded an understanding. We are going to demonstrate that it is 
not the case in the present situation. As we explain in the next section, the results of Golubev and Zaikin, being 
purely perturbative (1/t p is proportional to the first power of the fluctuation propagator), explicitly contradict to the 
straightforward calculation based on the perturbation theory. 



IV. MAIN IDEA OF THE CALCULATION AND RESULTS. 



The main idea of the calculation we ar e about to describe is the ob serv ation that in the strong magnetic field both 
results of the conventional t heory (2.45) and of the "new theory" (3.3) are purely perturbative in the interaction 
propagator U(Q,lu), see Eq. ( 2.25 ). Therefore, the truth may be revealed just by straightforward calculation of the 
interaction correction to the conductivity, see Fig. ^| We believe that we understand the qualitative physical picture. 
Nevertheless, we will not use this understanding in the calculation, trying, e.g., in the very beginning to identify 
"diagrams relevant for inverse inelastic time 1/Tj . . . " or to introduce cut-offs in the integrations using some ideas 
based on the qualitative understanding. Instead, we will take into account all the diagrams of the first order in the 
interaction and integrate over the whole interval of the energy transfer. The only small parameter we will afford to 
use is the inverse dimensionless conductance of the system. 
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FIG. 4. Diagrams for the first order corrections to the conductivity. Electron Green functions (solid lines) are exact for the 
given disordered potential. Wavy lines stand for the interaction propagator with Keldysh components C R , C A , C K , see Sec. |v| 

As usual, we start with expressing the correction to the conductivity through electronic and interaction propagators 
in a given disorder realization. This expression should be then averaged over disorder. For our purposes this averaging 
should be performed with taking the weak localization correction into account. It means that both zero and first order 
in the inverse dimensionless conductance 1/g should be kept while we average. Therefore, after the averaging and 
integration over the energies and momenta transfers, we will obtain corrections to the conductivity which are of the 
first and second order in 1 / g . 

It is important to emphasize that we evaluate the interaction correction to the observable (and thus gauge invariant) 
quantity - conductivity, rather than to a gauge noninvariant object like Cooperon. 

B efore we dig into details of the calculation, it is useful to realize what kind of correction we are interested in. In 
Sec. [IB the two kind of interaction effects - dephasing and interaction correction - were mentioned. Is it possible 
to separate the two corrections, when they are of the same order in 1/g? Fortunately, one can propose a procedure, 
which allows to perform this separation unambiguously. 

Indeed, let us write down the total correction to the conductivity in the form: 

5a = SawL + 8a c + <5c dcph + 5a C wh + ■ ■ ■ (4.1) 

where . . . stand for the corrections of the order higher than the second order in 1/g. Here, leading weak localization 
and interaction corrections to the conductivity are given by, see Eqs. (2.42) and Eqs. (2.28) 



5(Jwl = - 



In- 



1 



2ir 2 h tQ 



H 



5a c (T) 



2 h '" ( I t 



8a 



WL 



6*c(T) = - — 



fiD 



d = l, 



(4.2a) 
(4.2b) 



and time th S> h/T is given by Eqs. (|2 



The term Odcph appears as the first order correction to the time scale 
cutting off the weak localization correction. In accordance with the previous discussion, it should have the form 

2 

8a dc ph = —rr — ; d = 2; (4.3a) 



5a. 



dcph 



e 2 y/Dr H t h 



d= 1, 



(4.3b) 



where the scale 1/t v is to be determined. If the calculations of Golubev and Zaikin were correct, this scale would 
coincide with that given by Eqs. (77) and (81) of Ref. |. 

The last term in Eq. (4.1) characterizes the effect of weak localization on the interaction correction and vice versa, 
which do not break the time reversal invariance. On purely dimensional grounds they can be written as 



8a, 



CWL 



e 2 R a e 2 
8 a cwl 



a\ In Tt + ct2 In VLht + a.3 In Tr In fijjr) 



d= 1, 




(4.4a) 
(4.4b) 
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where oi is the conductivity per unit length of a one-dimensional system (wire), Rn is the sheet resistance of the 
two-dimensional film, and at are numerical coefficients to be found. One can separate contributions (^J^) and (4.4) 
using the apparent difference in the magnetic field dependence of these two terms. 

Now it is time to outline our strategy. We calculate the correction to the conductivity in the first order in fluctuation 
propagator. Next step is to average the result o ver d isord er, t aking weak localization into account in each term. Finally, 
we separate the result into terms of the form (|4.3|) and (4.4), thus evaluating the dephasing rate 1/t v . 

This program is described in details in Sec. |V[ but we would like to present and discuss the results right now: 



1. All the terms which scale according to Eq. (4.3) are produced by the expression of the form similar to Eq. ( 2.45 ) 
<W=^/ J ^/g) (^) ([C(Q,0)] 2 C(Q + g ,c)-|C(Q,c)| 2 C(Q + g ,0)), (4.5) 



where Cooperon C is defined by Eq. ( 2.45 ) and dimensionless function / is given by 



f( 



, , , de ( — tanh — | ( coth h tanh ■ 

V47V ./ \de 2T J \ 2T 2T 



(JtLY 1 

\2T) sinh 2 ^' 



(4.6) 



2. For the cross term of interaction and weak localization of the type (4.4) we found 



5a, 



CWL 



T 



n 2 h Tr 2 h<7d \KD 



1 

2d 



{2^Tr H ) 1 - d/2 ^j^r 2 (2 - d/2)C (2 - d/2) - ^( (3 - d) 



(4.7) 



where T(x) is the Gamma function and £(x) is the Riemann zeta-function. Here, we retained only most singular 
at d — > 2 contribution into the second not dependent on magnetic field term. 



Note, that in Eq. flf.Tj ) the ultraviolet convergence is provided by t he fu nction / (lo/T) from Eq. (4.6), which falls 
down exponentially at to ^> T (rather than by the limit ~ T as in Eq. (2.45)). The structure of / {oj/T) deserves some 
discussion. The coth term is the only term which was taken into account in Rcf. (7] This term indeed dominates for 
lu <C T. At high frequencies, however, this term is canceled out almost exactly by the second tanh termtil. It is this 
term, neglected in Ref. 0-^, that is responsible for the Pauli principle. The clearly indicates the inelastic origin in 
accord with the arguments of Sec. [n| The reason why the authors of Ref. missed this term is, as it is explained 
in Sec. VI, in the incorrect procedure of the disorder average. 

Since Eq. ( |4.5| ) is essentially the same as Eq. (2.45) in the frequency domain lu <C T, which gives the main 
contribution in low dimensional systems, this is not surprising that the results of the integration are perfectly consistent 
with Refs. p]j27|. In the leading order in Tth, we found 
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<5crdcph 



e 2 Tt h R a e 2 



2n 2 h h 

wh 



In ( T ™\ 
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which according to Eq. (4.3) means that 

1 TR a e 2 



h 2nh n \hJ 

1 T e 2 ^DTH~ 



d = 2: 



Ah hai 



d = l. 



(4.8a) 
(4.8b) 

(4.9a) 
(4.9b) 



Equations (4.9) are applicable provided that th <C t v ; in the opposite case 1/rjj in the right-hand side of these equa- 
tions, should be substituted by 1/t v and the procedure of Ref. ^7] is required for finding the temperature dependence 
for the magnetoresistance. 



The cross term of weak l ocal ization and interaction (4.7) is smaller than the dephasing term. We obtain for a one 
dimensional case from Eq. (4.7) 
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(4.10) 
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where o\ is the conductivity per unit length of one-dimensional system (wire), and C( x ) is the Riemann zeta- function, 
C(3/2) = 2.612.... 

Smallness of the term (R.10|) in comparison with the dephasing term (4.9) one more time illustrates the simple fact 



that the procedure of Ref. |27| is the parametrically justified program of summation of most divergent term in each 
order of perturbation theory in the interaction. 



The resulting conductivit y is given by the sum of Eq. (4.5) and Eq. (4.7). However, keeping the dephasing term in 
the leading approximation ( |4.S| ) and taking into account the c orre ction (4.7) would be overstep ping the accuracy of 
the approximation, since the contribution from oj ~ T in Eq. ( [4.5| ) is of the same order as Eq. (4.7). Therefore, the 
calculation of the integral in Eq. (4.5) should be performed with the same accuracy as Eq. (4.7). It gives 
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bumming Eq. @ and Eq. we obtain the final result for the overall effect of the electron-electron interaction 

on the weak localization correction 
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Taking the last term in brackets into account would be beyond the accuracy of the theory since the second order 
interaction correction gives the similar contribution. 



Let us write down the final results for the correction (4.12) in different dimensionalities 
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(4.13a) 
(4.13b) 



where C(l/2) = -1.461 . . ., C(3/2) = 2.612 . . .. 

To conclude, the results of the well-controlled perturbation theory contradict to those of Ref. 0. Thus, the statement 
made in this paper that the "analysis should in principle include all diagrams" is false. The theoretical conclusion 
about saturation of the dephasing rate at T — * contradicts not only the physical intuition but also the straightforward 
calculation. 



V. CALCULATION OF THE INTERACTION CORRECTION TO THE CONDUCTIVITY 



A. Conductivity for a given disorder realization 

In this section we will use units h = c = 1. We start with the_conventional expression for the current density j(r, t) 
in terms of the electron Green functions in Keldysh techniques 



j(r,t) = —- — Tr (fa-2fi)[V ri - V ra -2ieA(r)]G(ri,ti; r 2 ,t 2 ) 
4m 



ri = r 2 

t 1 =t 2 + = t 



(5.1) 
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where the trace is performed in the Keldysh space, the matrix Green function is defined as 



G(l,2) 



_{G R (1,2) G K (1,2) 



V 



G A (l,2) 



and the Pauli matrices in the Keldysh space are 











i 






i 



Entrees in Eq. (5.2) are given by 



2) = 9(h - t 2 )Wl)^(2) + ft(2)$(l)) 
iG A (l, 2) = -9(h - h)$(l)ft(2) + ^(2)^(1)) 
iG^(l,2) = <^(l)^t (2 ) _ ^t (2 )^(l)) 



(5.2) 



(5.3) 



(5.4a) 
(5.4b) 
(5.4c) 



where (. . .) means averaging over quantum mechanical states (do not confuse with the ensemble averaging !), 1 is the 
short hand notation for (ri,ii), ip(r,t) is the fermionic operator in the Heisenberg representation, and 9{x) is the 
step function. Spin indices are omitted and summation over them is implied whenever it is necessary. 
Let us now expand the Green function up to the first order in the interaction propagator 



£(1,2) 



C K {1,2) £ fl (l,2) 
£ A (1,2) 



(5.5) 



In the short-hand notation J d3 = J dt^dr^ the connection between Green functions in the presence (G) and in the 
absence (& n >) of the interaction can be written as 



6(1, 2) = «?<"> (1,2) + - f dZ fdA J2 G ( " ) (l,3)f l G (n) (3,4)f J G ( " ) (4,2)£ y (3,4), 



(5.6) 



For our purposes it will be sufficient to use the equilibrium values of the matrix elements of the interaction propagator 
£(1, 2). In this case these matrix elements do not depend on t\ + £2- Performing Fourier transform in the difference 
of the times = tj — t 2 



one obtains the following connection between the matrix elements^ 



C K (u) = coth — [L R (u) - L A (u>)] 



(5.7) 



In the leading order, the retarded and advanced parts of the interaction propagator depend on only on the difference 
of its coordinates. For the Coulomb interaction they are related to the screened interaction potential by 



C R (w,Q) = C A (-uj,Q) = -U(u>,Q) 



1 

27 



-iu + D(uj)Q 2 



(5.8) 



Here v is the density of states per one spin and U(u>,Q) is the dynamically screened interaction potential ( |2.26 ). 
To calculate the cross-term SacwL correctly, one has to include the weak localization correction to the diffusion 
coefficient D (see below). This is the reason why the diffusion coefficient D{u>) acquired the frequency dependence. 
The mesoscopic fluctuations of the propagator C(u>) are disregarded in the present calculation, since they lead to the 
higher order in 1/g correct ions to the co ndu ctivity. 

Now, we substitute Eq. (5.6) into Eq. (5.1). In the equilibrium, the longitudinal current vanishes, and therefore all 
the Green functions should be expanded up to the first order in the external vector potential A cxt which represents 
the applied electric field E = d t A ext : 



G (n) (l,2) = G (0) (1,2)+ / d3 G (0) (l,3) jA cxt (t 3 ) G (0) (3,2). 



(5.9) 



17 



The current operator in the coordinate representation is defined as 



fi{r)jh(r) = [(V/0 h /iV/ a ] - ^/jW/^r) 



(5.10) 



for any two functions fx %. 



The time Fourier transform of the matrix elements of equilibrium Green functions G^ in Eq. (5.9) are given by 



G K (e; rx,r 2 ) = tanh JL [G R (e; n,r 2 ) - G A (e; r ls r a ) 



G« (e;n,r a ) = 
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e + /x - if + iO 



G- 4 (e;ri,r a ) = 



H-H-iO 



(5.11a) 
(5.11b) 



where H is the exact one-electron Hamiltonian for a given system, and \x is the chemical potential. In Eq. (5.11) 
and in all subsequent formulas we will omit superscript (0) since only such equilibrium Green functions appear in all 
further calculations. 

Substituting Eqs. ( |5.9| ) into Eq. Q5.6D and the result into Eq. (5.1) we obtain after simple algebra the Ohm law 
j = aE with the conductivity tensor 
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(5.12c) 
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where V is the volume of the system, j a is the compone nt of cu rre nt ope rator defined in Eq. ( [5.1C ) , and the short-hand 
notation Gy (e) = G(e;r,,rj) is introduced. In Eqs. (5.12b) - (5.12c) only the symmetric part of the conductivity 
tensor is retained. We emphasize that all these formulas are exact in the first order of perturbation theory, i.e. 
we made no assumptions about the behavior of the Green functions or the rele vant f reque ncies in the interaction 
propagator. At this stage separation of the interaction correction into two pieces (|5.12c| ) and ( pM2d| ) is just a matter 
of convenienc e phy sical meaning of which will become clear shortly. 

Equations ( 5.12 ) , which express the conductivity in terms of the exact retarded and advance Green functions for a 
given realization of potential, are the main results of this subsection. In order to find the physical conductivity of a 
large sample, one has to perform the disorder averaging. This is the subject of the following subsection. 



B. Averaging over disorder 



Let us start the disorder averaging of Eqs. ( 5.121: ) - ( 5.12d ) with the observation that any mean product of the 
electron Green fu nctions does not depend on their common energy e. This fact allows to perform the energy integration 
in Eqs. ( 5.12b| ) - ( 5.12d ). It is also important and useful to notice that after the averaging, t he p roducts containing 
only retarded or only advanced Green functions vanish. Using this feature, we simplify Eqs. ( 5.1 2j ) to the form 



'a/3 



1 

4tt 



drxdro, 



((G 



L2 ( e ) 3 aG 2\ 



(5.13a) 



18 



deph 

a/3 — 



I 



dridr 2 dr 3 dr4 
V 




2Tsinh^ 



[C& (to) - (w)] x 



2T , 



{<(7«Gf a (e) j>G& (e) (e - w ) (e))) + (e) j>G& (e) Gf 4 (e - w ) G« (e))) 

+ ((j tt Gf 3 (e) G* (e - to) j G$ 4 (e - co) G* (e))) + 



(5.13b) 



1 /" dridr 2 dr 3 dr4 
Im / 



2tt 



(lj coth — ^ 



2TV 



{ 



(O a Gf 2 (e)i>G« ( £ ) G3I ( e - W ) G« (e)» + ( (ja G12 (e)j Gi 3 (e) G* 4 (e - w ) G 41 (e)» 
-«G& (e) j a G& (e) (e) G 4 « (e - u;)» - «G& (e) j Q G^ 4 (e) G% (e - W )^G& (e - u))) 

+2((Gf 1 (e)3 Q G£ (e) G% (e - u) jpG^z (e - w))) + a ~ /?} 



(5.13c) 



where ((...)) stand for the disorder average. 

Equations ( 5.13 ) express the average conductance in terms of the average of products of Green functions of a non- 
interacting system. Those averages should be performed with_taking the weak localization correction into account. 
We will perform the average using the diagrammatic techniquetlH; the same results can be obtained with the help of 
supersymmetric non-linear cr-modelo. 



Evaluation of the mean conductivity of non-interacting electrons a m , i.e., averaging of the equation (5.13a), is now 
a textbook problem. Summation of the diagrams in Fig. |5| gives a a /3 = S a fi (ff + 5^wl) where Drude conductivity a 
and the weak localization correction given by (1.1) and the weak localization correction 6o~wl are given by (1.1) and 
(2.5) respectively. 
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FIG. 5. Diagrams for the averaging of a nt . (a) Drude conductivity, (b) Weak localization correctionQ; (c) More convenient 
for the further usage representation of diagram (b). Hereafter, solid lines denotes disorder average Green functions and dashed 
lines denote disorder. 

Before averaging the interaction correction, it is convenient to Fourier transform Eqs. (5.13b) and ( [5.13c ) over the 
difference r 3 — r 4 . We obtain from Eq. ( 5.13b| ) 
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where hereafter we use the short hand notation 
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for any functions Fw. 

Three last terms in brackets of Eq. ( 5.13c ) can be simplified with the help of the identities 

J dr 3 G R (e, n, r 3 ) j a G R (e, r 3 , r 2 ) = -ie (n - r 2 ) a G R (e, n,rj) ; 
dr 3 dr 4 G R (e, n, r 3 ) j a G R (e, r 3 , r 4 ) j^G^ (e, r 4 , r 2 ) j> = 

- ~ (n - r 2 ) Q (r x - r 2 ), G* (e, n, r 2 ) + ^^G* (e, nra) , 

which follow from the gauge invariance of the theory and do not require any translational invariance, see Appendix. 
We find 
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Deriving Eq. ( p.lq ) we used the fact that 
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Now we are in the position to perform actual disorder average of the products entering into Eqs. (5.14) and (5.16). 
Diagrams for all the terms from Eq. ( 5.14 ) are shown in Fig. ^|. They can be obtained by cutting the corresponding 
Green functions on the diagrams of Fig. |b. The analytic expression for these diagrams is 
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where the Cooperon pr opa gator is defined by Eq. ( 2.45| ). Substitution of Eqs. (5.17) and ( j5~8| ) into Eq. (5.14) 
immediately yields Eq. (4.5). 






FIG. 6. Diagrams for the averaging of products (a) (GcG&jpG&G&G&fiui (b) ((jaG^jpG^G^G^)):^, and (c) 

( (ja Gf 3 G 3 2 j 13 G 24 G 41 )) 34 

Diagrams for the averaging of the terms entering into Eq. (5. IE) are a little bit more complicated, see Figs, fj] - 1|. 
The analytic expression for the diagrams in Fig. [7] is 
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where Cooperon propagator C{— u>, Q) is defined by Eq. ( p.45| ) and the diffuson is given by 
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-iuj + DQ 2 



(5.18) 



(5.19) 



The first term in Eq. (5.18) is the leading classical approximation, Fig. ^a and all the other terms come from the 
interference corrections which are shown term by term in Fig. pp-d. 
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FIG. 7. Diagrams for the averaging of products {{j a Gi 2 jpG 2 :iG 3 4G 41 ))34 



The third term in brackets in Eq. (5.1£) is the density - density correlation function with the well known answer, 
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see Fig. (|) 
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where the last term is nothing but the weak- localization correction to the diffusion constant. 
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FIG. 8. Diagrams for the averaging of density-density correlation function ({Gi|G2i))i2. 



Finally, the last term in brackets in Eq. ( 5.1£| ) is averaged according to Fig. ^. The analytic expression for this 
diagram is 
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FIG. 9. Diagrams for the averaging of product {(G^jaG^G^)) <23 



Now, we are ready to calculate the inter action corre ction (5.16). First we take i nto ac coun t leading 1/g approxima- 
tion only, i.e. retain t he fi rst terms in Eq. ( |5.18 ) and ( 5.20| ), neglect terms ( 5.17a ) and fl5.21 ), and use the i ntera ction 
propagator C A from (5.8) with the bare diffusion constant D. It immediately gives the known result, Eq. (2.27). 

In order to find the cross-term of weak interaction and weak localization, we have to take into account not only 
all the remaining terms in Eqs. ( |5.17a ), (5. IS), ( 5. 20| ) and (5.21) but also the weak localization correction to the 
interaction propagator (|5.8[) in the leading approximation 
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(5.22) 



where the last term i s the w e ak lo calizat ion c orr ection to the dif fusio n constant. 

Substituting Eqs. ( 5.17a ), ( [5.18 ) and ( 5.20 ) - ( 5.22 ) into Eq. ( 5.16| ) and keeping only terms of the second order in 
1/g, i.e. only those involving two integrations over the wavevectors Q, q, we obtain 
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Derivation of Eq. (4.7) from Eq. (5.23) is now just a matter of a lengthy albeit straightforward calculation. We will 
write most important intermediate formulas below. Integration over the wavevectors is performed first. To carry out 
this integration, it is convenient to represent each Cooperon and diffuson propagator by additional time integral 
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After this manipulation, the integrals over wavevectors become Gaussian and can be immediately calculated. We find 
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(5.25a) 



where r(x) is the Gamma function. Equation (5.25a) is written with an accuracy up to term of the order of 

uj d ~ 3 (ujth ) d ^ 2 ■ This accuracy however is sufficient for our purposes, since all the Cooperons and diffusons here 
are functions analytic in the complex semiplane Imw < 0, and therefore integration over lo is determined by u > T. 
On the other hand, we are considering the case Tth ^> 1, and all the corrections to Eq. (5.25) are smaller by this 
parameter. Other integrals calculated in the same approximation are given by 
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Substituting Eqs. (5.25) into Eq. ( 5.23| ), retaining only most singular at d — > 2 coefficients in front of u d 3 , and 
omitting indices a,/3, we find 



2irv 2 (± n D) 
Using formula 



rim 



duj 
2^ 



d 

du> 



Im 



(u) coth 



2TJ 



4r 2 (2-d/2) . d/2 _ 2 / 1 ^Z 2 " 1 

j M 1 - 



4 , d -3 



2-d 



(5.26) 



dx d 
(ia;) Q da; 



fx coth — 
V 2 



(2tt 



2a , s 
— C(a), 



with being Riemann zeta-function, we obtain Eq. (4.7). 



VI. DISCUSSION OF THE RESULTS 



A. Theory 



In Sec. we demonstrated explicitly that dephasing can be caused only by e-e collisions with the energy transfer 
smaller or of the order of temperature, see Eq. (45). This conclusion clearly contradicts the results of Ref. ||. It 
means that the procedure proposed in Ref. fails already on the level of the first order perturbation theory in the e-e 
interactions. This calculation is sufficient to make the conclusion that the "old" rather than "new" theory describes 
the temperature dependence of the dephasing rate correctly and thus to end the pointless discussion. However, the 
errors in Ref. ^| originate from a quite typical misuse of the semiclassical approximations. For this reason we would 
like to highlight thes e erro rs and to di scus s them in more detail. 

According to Eqs. (2.45), (4.5) and ( |3.3| ), dephasing can be found by determining the coefficient in front of the term 

2 

[C(Q, 0)] . The Q — ► O and r — > oo limit of this coefficient has a meaning of the inverse dephasing time while 
its Q and independences represent renormalization of the diffusion constant D and of the residue of the Coo per pole 
respectively. These renormalizations contribute eventually to the in terac tion correction to the conductivity (^?]) . 
Therefore we can focus only on the first term in brackets in Eq. ( 5.23 ): 



o f duj 
-1m / — 

2ttv 2 / 2tt 



d ( uj 
— uj coth —— 
duo \ 2T 



d d qd d Q 



(27T 



, d C a {uj, q)C 2 (0, Q) [C(-w, Q + q) V(-u, q) 



(6.1) 



Note that in the corresponding equation of Ref. |9|, see Eq. (3.3), only the first Cooperon term in the brackets is 
present. The second diffuson term changes the result dramatically since at Q — > O and Tp — ► oo it simply cancels the 
first one. 

In order to understand why the diffuson term and thus the cancellation are missed in Ref. [|, let us examine step 
by step the derivation proposed in that paper. 

1. Authors of Ref. [| used Keldysh technique and decoupled the interaction by Hubbard-Stratonovich transfor- 
mation. This step is rather standard and the result in the perturbation theory should be consistent with our 
Eq. E/B. 



2. Their next step was to write down the formally exact path integral expression for Green functions in a disordered 
potential 



G R (e ]ril r 2 ) = 



G A (e;r l7 r 2 ) 



die 1 



r(t)=r 2 

r(o)=7"i 
r(t)=r 2 

r(o)=ri 



V[r(t)]exp l-S[r(t)} 



Z>[r(t)]exp(--S [r(t)) 



(6.2) 



where S [r(t)) is the classical action over the trajectory r(t). 



25 



3. Afterwards an averaging of products of Green functions was performed. During this average only some paths 
were selected by hand. It is this selection that, as far as we understand, causes all the errors. 

Indeed, consider the average 

({G R {e;r 1 ,r 2 )G A (e;r 1 ,r 2 )))=2^V(0;r 1 ,r 2 ), n, r 2 ) = J ^ (6-3) 

Here, for simplicity, we consider the case of zero magnetic field and do not distinguish between diffusons and Cooperons. 
This average can be obtained semiclassically by pairing the paths shown in Fig. |lO| a and it is reproduced correctly in 
Ref. ^. The same is true for the simple average involving four coordinates, see Fig. |To|b: 

{{G A {e- n, r 2 )G H (e; r x , r 3 )G R (e - u; r 3 , r 4 )G /? '(e; r 4 , r 2 )» = -2nuV(0; n, r 3 )V(-Lj; r 3 , r 4 )P(0; r 4 , r 2 ). (6.4) 

This expression corresponds to the diagram @a. 

However, the other higher order products involve intersection of the classical paths, see Fig. [Io| c. Such configurations 
contribute in, e.g., average of the type 

((G R (e; ri ,r 2 )G A (e; r 1} r 3 )G R (e - u; r 3 , r 4 )G A (e; r 4 , r 2 ))) (6.5) 

= 2ttv f dr 5 {iuV(0;ri,r 5 )V{-u]r 5 ,r 3 )V{-Lj]r 5 ,r4,)T>(0-,r2,r 5 ) 

+ D [V 5 2?(0; n, r 6 )V{-u; r 5 , r 4 )] [V B 2?(-w; r 5 , r 3 )V(Q; r 2 , r 5 )] 
+ 2D [V 5 P(0; n, r 5 )\ [V 5 V(-u; r 5 , r 4 )] 23(-w; r 5 , r 3 )P(0; r 2 , r 5 ) 
+ 2D D(0; n, r B )X>(-w; r 5 , r 4 ) [V B 2?(-w; r 5 , r 3 )] [V 5 P(0; r 2) r 5 )]} , 

see diagram in Fig. fjjp. The point r 5 has the meaning of the intersection point between classical trajectories. These 
intersections were not taken into account at all in Ref. [| 

After performing the Fourier transforms over the coordinate differences r 3 — r 4 and ri — r 2 , one finds 

& - 



[DQ 2 ] 2 iuj + D(Q~q) 2 
j— i 2ttv 2-kv ini/ DQq 

(p-5D 



[DQ 2 ] 2 [iuj + Dq 2 } [DQ 2 }[ioj + Dq 2 } 2 [DQ 2 ] 2 [iuo + Dq 2 } 2 

where the last term will vanish after the angular integrat ion. One can see that signs of Eqs. ( |6.3| ) and ( |6.4| ) are 
opposite, which eventually leads to the cancellation in Eq. (|6.1|). 
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FIG. 10. Example of the classical paths contributing to the averages of products of the Green functions. 

In fact, taking into account the switches between intersection of the classical paths is very crucial: mistreatment of 
such switches leads, e.g. to the violation of the conservation of number of particles and all the theoretical predictions 
of similar quality! In the diagrammatic technique such switches are described by so-called Hikami boxes (the shaded 
boxes in Fig. |t] - ||). Semiclassical approach usually has difficulties in the description of the transitions between 
the classical trajectories. Physical reasons for this difficulty were first discussed by Larkin and OvchinnikovEj in 
application to the superconductivity, and in Ref. [36] in application to the weak localization. 

B. Experiment 

Let us compare the available experimental results with predictions of the theory by Golubev and Zaikini. Authors 
concluded that "according to Eqs. (81), [of Ref. |j| in 2d and 3d systems the decoherence time becomes independent 
on T already at relatively high temperatures. In the 3d case such temperatures can be of the order of the inverse 
transport time, while for 2d systems we haveT, - (2rln(p|aZ)) ".If we now recall that the momentum relaxation 
time t in the weak localization experiments with disordered metals is typically 1CP 16 — 10~ 14 sec, we will have 
to believe that the phase-breaking time due to electron-electron scattering should be T-independcnt in 2d and 3d 
conductors at all temperatures below the melting point! This estimate for t v is also much smaller than the electron- 
phonon scattering time r e _ p ^ in the temperature range T < 20X, which is typical for the experimental studies of 
the quantum corrections to the conductivity. In this situation, the T-dependent weak localization corrections to the 
resistance of 2d and 3d conductors should be non- observable. For Id and 2d conductors, the strong-localization regime 
should be non- observable either. Both corollaries of Eqs. (81) of Ref. ^ are in a profound contradiction with the bulk 
of experimental data, which we review below. 

1. Three-dimensional conductors. 

The study of the WL corrections to the conductivity of 3d systems has been focused mostlyon thick (a ^> L v ) 
disordered metal films, metal glasses, and heavily doped semiconductors (for reviews, see Refs. p7H39[). For all these 
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systems, the well-pronounced temperature dependence of the WL corrections has been observed. From the analysis of 
the WL magnetoresistance, the temperature dependence of, the phase-tweaking time has been extracted. The typical 
T<p(T) dependences, obtained for 3d disordered films of Cuc3 and Cu-Gecil, are shown in Fig. [ll]. The electron-phonon 
scattering governs the phase-relaxation in these 3d conductors over the whole temperature range T = 1.5 — 4QK. The 
solid line T e _ p ^ cx T~ 3 corresponds to the electron-phonon relaxation time due to the longitudinal phonons in pure 
Cu (see, e.g., Table 3 from Ref. |37]). The agreement with the experimental data can be furth er .| im proved if one takes 
into account the contribution of transverse phonons, which dominates in disordered metalsclle3_ The temperature 
dependence of r v (3d) due to the electron-electron interactions with the small momentum transferd 



V2 (k B T\ 3/2 

(y is the density of states) is shown in Fig. [ll] with the dashed lines. The upper limit of t v accor ding to Ref. | 

h 



(6.6) 



(6.7) 

(<7 is the conductivity, and r is the elastic scattering time) is shown in Fig. [ll] with the dotted lines; this estimate is 
much smaller than the experimentally measured time. The discrepancy here is as large as 5 (!) orders of magnitude. 
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FIG. 11. The temperature dependence of t v in 3d metal films: O ~ disordered Cu films (1) with p = 6.6 ■ 10 

s, Ref. 0; A — Cu . 9 Ge .i films (2) with p = 2.8 • l(T 5 fi ■ cm, D = 9.3 cm 2 /s, r = 1.1 



L> = 3.9 cm 2 /s, r = 4.8 ■ 10' 
Ref. fij]. The solid line: r e - ph 
respectively. The dotted lines: the maximum t v (3d) according to Eq. (6.7) 



5 n 
■ 10 



cm, 



(pure Cu) = 2.2 • 10~ 7 /T d s. The dashed lines: r ee (3d), Eq. (6.6) for these samples 1 



and 2 



Similar discrepancy between the experimental values of t v and predictions of Golubcv and Zaikini has been observed 
for metal glasses (see, e.g. Refs. P8| , [43| -^5|). For example, for {Ca$oAl2erk)7.iAu2.9 with p = 1.3 • 10~ 4 Q ■ cm and D ~ 1 

ne could 
reported 



cm 



/s the experimental value of t v (2.5K) = 7.5 • 10 10 s was observedEI On the other hand, from Eq. (6 



expect t v 
TJ4.2K) 



< 4 
~ 2 



10 
10 



-14, 



s. For Y%oSi2o with p = 5 • 10 
11 s, whereas Eq. (6.7) gives t v c 



Q ■ cm, t 
3 • 10- 16 s. 



1 ■ 10~ 16 s and D ~ 0.7 cm 2 /s, Bieri et al. 



2. Two-dimensional conductors. 



Ti 
film 



Tie. weak localization regime. — Figures [l2] and, [l3| show the typical temperature dependences of t v in 2d metal 
sc!TC3 and GaAs — AlGaAs heterostructuresEll, correspondingly. At low temperatures (T < 1 — 5K), the phase- 



breaking is governed by the quasi-elastic electron-electron scattering (the Nyq uist p hase breaking) . The experimental 
values of t v are usually half as large as the estimate of T 9 (2d), given by Eq. ( 2.381 ). 

On the other hand, the experimental data for t v exceed by far the upper limit of T v (2d) according to Eq. (81) of 
Ref. 



t gz {2D) 47rt 2 t' 



(6.8) 



The discrepancy of Eq. (3.8) with the experiment is especially 
of magnitude for cryo-deposited Mg films at T = OAK, Ref. [ 



sronounced for the most disordered samples (4 orders 
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FIG. 12. 
r = 4.6- 10 

2 



Tbei temperature dependence of t v in 2d metal films: 1 (Q) - Mg filra-jwith Ra = 22.3 Q, D — 3.8 cm 2 /s, 



2 (AjU- AZ film with Ra = 112 fl, D — 10.4 cm 2 /s, r = 1.7 • 10~ 15 s; El 3 (□ ) - Si film, R a = 630 fi, D = 8 



cm 2 /s, r = 2.8 ■ 10 14 s.c3 The solid lines: the t emp erature dependences of r v {2d) [Eq. (2.38b)] for these samples. The dashed 
lines: the maximum Tcz(2d) according to Eq. (6.8). 



The strong localization regime. With increasing the sheet resistance of 2c? samples up to ~ 30 k£l, the crossover to 
the strong localization regime has been observed for both quasi-2rf metal films and 2c? electron gas in semiconductor 
structures (see, f.L, Refs. |50j-|54|). A little extension of the Golubev-Zaikin logic leads one to the conclusion that the 
"zero point motion" also suppresses the interaction correction to the conductivity since the latter is also caused by 



the interference, see Sec. [IB. As the result, the conductivity should be temperature-independent at T < Tqz — 
fo/ T Gz(2d). For all the samples studied, Tqz is very large; for example, the samples of Refs. pl| -|53| are characterized 
with Tqz — 10 5 i^. On the other hand, all the crossovers to the strong localization were observed at temperatures 
of the order of IK — 10K. These experimental data rule out the statement by Golubev and Zaikin that ". . . strong 
localization does not take place at all and the lc? and 2c? metals do not become insulators even at T ~ 0" . 
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FIG. 13. The temperature dependence of r v in 2d GaAs — Al x Ga\ 



cAs heterostructuresl 

-2 



E 



1 (□) 



cm 

-13. 



i = 0.87 ■ 10 1 

R D = 4.3 kn, D — 51 cm 2 js, r = 6.3 • 10~ 13 s ; 2 (O) - n = 2.86 ■ lO 11 ^" 2 ; i? D = 4.0 fcfi, £> = 55 cm 2 /s, r = 2.1 • 10 _ "s; 3 
(A) — n = 7. 11 ■ 10 11 cm~ 2 ; = 300 fi, D — 753 cm 2 /s, t = 1.1 • 10~ 12 s. The solid lines: the t emp erature dependences of 
rjv(2D) [Eq. ([2.38b)] for these samples. The dashed lines: the maximum T v (2d) according to Eq. (6.8). 



3. One-dimensional conductors. 



In Ref. |TT|. Golubev and Zaikin compared their results with the experiment by Khavin et aL@ and claimed a 
convincing agreement. We will show below that_in fact, the agreement exists only in the temperature region where 
their result coincides with the conventional oneEH. Outside this region the comparison becomes meaningless because 
of the observed crossover to the strong localization regime. nn 

The weak localization regime. — In one dimension, the expression for T v (ld) obtained by Golubev and Zaikinolj, see 
Eq. (77) of Ref. | or Eq. (10) of Ref. §, is 



h e 2 2D 

TCzO-d) 7TO~l V T 



2k B T ; 

1 H r — a/tgzT 



(6.9) 



In the temperature range fcflT> %/ \Fhpf, this expression coincides up to a numerical factor of the order of unity 
with the Nyquist phase-breaking time in ldj23 see Eq. (2.38a). 

It has been established experimentally that the phase breaking in Id conductors (both metal wires and semiconduc- 
tor structures) is well described by T v {ld) over a wide temperature range (typically, below several Kelvin) (see, e.g. 
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Refs. ^5|-|60|). Thus, no wonder that the expression (6.9) is in a ccord with the experimental data in the regime where 
the "new theory" coincides with the old-fashioned Eq. ( |2.38a ). We would like to point out again that the explicit 
comparison of 1/tqz with the experiment is meaningful only if the analytic expression for the m agnetoresistance is 
provided. Since Golubev and Zaikin were not able to obtain such formu la, similar to Eq. ( 2.42c ), the procedure of 
comparison of their results with t v extracted with the help of Eq. (2.42c) is at least problematic. 

On the other hand, in the temperature range ksT <C Ti^/t^t, Golubev and Zaikin predict saturation of the phase 
breaking length: 



L 



GZ 



3 
2^2 



1/4 



1 



N 



(6.10) 



where N is the number of the transverse channels in the wire, Lqz — \/Dtgz, and £ is the Id localization length. 
Once again, we do not know what the numerical factor in th is f ormula means, since Golubev and Zaikin did not 
derive the explicit formula for the conductance. [In fact, Eq. ([?!]) seems to be erroneous even within the Golubev- 
Zaikin logic: the contribution to the zero-temperature term comes from the wave-vectors of the order of the inverse 
clastic mean free path 1/7. Therefore, to be consistent with the "new theory", one "should" use for a wire with the 
cross -sec tional dimensions greater than the elastic mean free path either two- or three- dimensional formulas (|6.8| ) 
and (6.7), depending on the structure of t he w ire.) 

Anyway, the author s of Re f. |ll"| use Eq. (6J3) for the a nalys is of the data. In the regime where deviations from the 



conventional formula (2.38a) can be expected from Eq. (6.9), the information on t v can not and was not obtained 



from the experiment, because the samples were already in the strong localization regime. 

The strong localization regime - It has been shown recently that the Id wires fabricated from .Si 5 -doped GaAs, 
demonstrate the crossover from weak to strong localization with decreasing the temperatureEiJ£3. This crossover is 
driven by both localization and interaction effectst3; in particular, the importance of the localization effects has been 
demonstrated by observation of (1) the shift of the crossover by the weak magnetic field; (2) the orbital magnetoresis- 
tance on the insulating side of the crossover, and (3) the doubling of the localization length in strong magnetic fields. 
These experimental facts argue against the conclusion of Golubev and Zaikin "strong localization does not take place 
at all and the Id and 2d metals do not become insulators even at T — 0" . 



VII. CONCLUSION 



This paper is devoted to the quantum transport of electrons in disordered conductors in the weak localization 
regime. The motivation to write a paper on the subject that was developed almost 15 years ago was caused by 
recent experimental and theoretical attempts to reconsider the existing theory of destruction of the one-particle 
phase coherence (dephasing) in a system of degenerate interacting fermions. It was proposed that, in contrast to basic 
conclusions of this theory, the dephasing rate remains finite even at zero temperature due to the zero-point oscillations 
of the electric field. 

We criticized this proposal in details from several points of view. 

1. We recalled and discussed in details the physical picture of the quantum corrections to the properties of the 
disordered conductors and tried to demonstrate that such a revision looks absurd from the qualitative point of 
view. 

2. We made a straightforward calculation of the first order in the e-e interaction correction to the conductivity 
at low dimensions in stron g magnetic fields, taking into account all the relevant diagrams. Results of this 



calculation, see Eqs. ( 4.13 ), are in a perfect agreement with the "old" theory and explicitly contradict the 



zero-point dephasing calculations. 

3. We highlighted the sources of the mistakes of the new "theory" . The errors were found to be due to the misuse 
of the semiclassical approximation in disorder averaging. 

4. We reviewed existing experiments on magnetoresistance in 3d, 2d, and Id dimensions, and demonstrated that 
the measured dephasing rate is sometimes more than five orders of magnitude smaller than the limit, which 
follows from the zero-point-motion "theory" . 

We really hope that arguments presented are sufficiently strong to stop the "theoretical" discussion of the zero 
temperature dephasing. 

Our paper also contains the new explicit formulas for the effect of the interaction on the weak localization due to the 
clastic scattering. Such effects, are shown to be 1/g 2 corrections which are smaller than the main terms considered 
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by the conventional theory everywhere except the vicinity of the met al-insulator crossover. We should warn the 
practitioners of the field that in addition to the extra terms in Eq. ( 4.13 ) the first order correction to the conductivity 
in the Cooper channeB and the second order weak localization correction to the conductivity^ has to be taken into 
account for comparison with the experiment. 
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APPENDIX A: DERIVATION OF IDENTITIES (5.15) 



Let us perform the gauge transformation of the Green function 

G R (e, ri,r 2 ) -► e ^Mri-r 2 ) G R ^ ^ Ti 



(Al) 



where A is an arbitrary vector independent of the coordinates. Under this transformation, the exact one-electron 
Hamiltonian changes according to the rules 



H -> H -jA 



s 2 A 2 
2m 



(A2) 



where j is the current operator ( 5.1 Cl| ) . 

Therefore, it follows from the definition of the Green function that the relation 



eA-(r 1 -r 2 ) 



1 




1 


.e-H + i0_ 


12 


_e-H+jA- 



e^A 2 



iO 



(A3) 



holds. Expanding both sides of Eq. ( |A3| ) up to the second order in A, we obtain Eqs. ) 
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